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Abstract I review recent results concerning the shape of drifting subpulse pat- 
^ ' terns, and the relationship to model predictions. While a variety of theoretical 

CO , models exist for drifting subpulses, observers typically think in terms of a spatio- 

temporal model of circulating beamlets. Assuming the model is correct, geometric 
parameters have been inferred and animated "maps" of the beam have been made. 
f^ ' However, the model makes very specific predictions about the curvature of the drift 

VO i bands that have remained largely untested. Work so far in this area indicates that 

^P ' drift bands tend not to follow the prediction, and in some cases discontinuities are 

r"| . seen that are suggestive of the superposition of out of phase drift patterns. Recent 

Q^' polarimetric observations also show that the drift patterns in the two orthogonal 

polarisation modes are offset in phase. In one case the pattern in one of the modes 
^ ' shows a discontinuity suggesting no less than three superposed, out-of-phase drift 

C/3 , patterns! I advise caution in the interpretation of observational data in the context 

^ ' of overly simplistic models. 
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1 THE CAROUSEL MODEL 



Although the drifting subpulses phen omenon has now been approached from a number of 
different theoretical perspectives (e. g. Zheleznvakov 1971; Kazbegi_et_aL[Jl99]l IWrighlJ l2003t 
IClemens fc RosenI 120041: lYoundl2004£ .Goeoberidze et al...2005: .Fung et alJl2006|) . most obser- 
vatio nal studies have assumed the applicability of the so-called carousel model iJRudermanl . 
Il972|) . This model postulates that the modulations have a spatio-temporal origin: the subpulse 
structure is due to the passing of the line of sight over multiple "beamlets" , while the temporal, 
pulse-to-pulse modulation is due to the slow circulation of the beamlet system (Fig. ^^a)). In 
the usual physical model, the beamlets consist of radiation beamed tangentially to local mag- 
netic field lines, within "tubes" of plasma, each fio wing outward from a localised breakdown 
( "spark" ) of a potential gap over the magnetic pole l|RudermanL Il972|) . The circulation of the 
sparks is attributed to E X B drift. 
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Fig. 1 Diagram of sub-pulse emission geom.etry, reproduced from 

lEdwards fc Stapoersl l)2002l) . In part (a) the angular momentum, magnetic mo- 
ment, and line-of-sight vectors are shown with symbols L, /x, and d respectively. Part 
(b) shows the main vectors as they appear in the plane they share when the sight-line 
makes its closest approach to the magnetic pole. The angles C, (between the sight-line 
and the spin axis), a (between the magnetic and spin axes) and /3 = ^ — a (between 
the sight-line and the magnetic axis at their closest approach) are shown. In part (a) 
a second line of sight (d') is drawn, along with the meridian it shares with with the 
magnetic pole. The spherical triangle formed by L, fi and d' is duplicated in part (c). 



Edwards &: StappersI l|2002(l showed that the carousel model makes very specific predictions 
regarding the shape of drift bands. Specifically, assuming that the beamlets are uniformly spaced 
in magnetic azimuth, the variation of subpulse phase with pulse longitude obeys a geometric 
relationship. This is because it is linearly tied via the number of beamlets (N) to the magnetic 
azimuth of the observer (ip), which in turn obeys the following relation: 



tantp 



sin (J3 sin ^ 



cos ^ sin a — cos sin (^ cos a ' 



(1) 



where 4> is pulse longitude, ^ is the angle between the spin axis and the line of sight, and a is 
the angle between the magnetic and spin axes (see Fig. ^b)). The reader may recognise the 
similarity between this relation and the common "rotating vector model" for the position angle 
of linear polarisation, x- The reason for this similarity is that all of the pertinent angles are 
related according to the spherical triangle of Fig.^c). The subpulse phase, 9, is linearly related 
to magnetic azimuth, plus a generally small additional term describing the circulation of the 
carousel over the course of the pulse: 
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Here 0' is an arbitrary constant, Pi is the pulse period, P3 is the observed pulse-to-pulse 
periodicity of the drift bands (including a sign which specifies the direction of the drift slope), 
and n is an integer specifying the order of aliasing present. 

It is important to note that by examining the dependence of modulation phase upon pulse 
longitude, this method tracks the drift band peaks according to local maxima along lines of 
constant longitude. This is in contrast to the more common technique of finding local maxima in 
individual pulses (i.e. lines of constant pulse number). The former method offers an important 
advantage: the amplitude-phase decomposition conveniently separates beamlet azimuthal spac- 
ing (i.e. circulation or subpulse phase) from colatitudinal amplitude windowing (i.e. subpulse 
amplitude and mean profile shape). In the latter method, subpulse peaks are "pushed" in pulse 
longitude, in the direction of increasing s ubpulse amplitude, and the refor e do not cleanly track 
the azimuthal spacing of beamlets (see Edwards & Stapperal bood and Edwards fc StapDera 
EqOS;). 



2 WHY TEST THE CAROUSEL MODEL? 

While the carousel model has been in use for over two decades in the interpretation of drifting 
subpulse patterns, until recentl y very little u se had been made of the opportunity presented 
to test the model (however, see lWrightlll98l|) . Many examples exist of inferences being made 
on the basis of the correctness of the model, without any checks being made of the validity 
of this assumption. For brevity I will focus upon a critical look at one particularly interest- 
ing seri es of studie s, namely the industry of "polar ca p mapping" of syste ms of drif ting sub- 
pulses i DesIroande fc^ankii^, ^1999,; ^Deslrom Rankin . 2001: Asgekar fc Deshpande. 2001i 

iRankin et al.U2003HAsgekar fc Deshpandd. 12005(1 



lOeshpande fc RankinI (|l999j (hereafter DR99) detected, for the first time, a pair of "out- 
rider" components surrounding the first harmonic of the drifting subpulse pattern in the fluc- 
tuation spectrum of PSR B0943-K10. While relatively marginal in terms of signal to noise ratio, 
and only detected in a small section of one of several observations, the detection was never- 
theless exciting since it was potentially supporting evidence for the carousel model. This is 
because under that model, any overall overall pattern of intensity around the ring of beamlets 
would result in the convolution of all Fourier components with the Fourier transform of the 
azimuthal intensity dependence ([Edwards & StappersH2002l) . For an intensity pattern consist- 



ing of a strong steady component plus a weak variation oc sin(^ — V'o), this adds a small pair 
of components spaced I/P3 from the parent component. Although no such "outriders" appear 
around the DC component, DR99 pursued the possibility that this was the explanation for the 
extra components around the first harmonic. 

Having assumed the applicability of the carousel model, DR99 set out to establish the values 
of all geometric parameters necessary to invert the transformation between the magnetic frame 
and the observer's frame, in order to "map" the observed intensity onto an animated display 
of the beam intensity distribution (or, equivalently, the polar cap excitation distribution). The 
first hurdle in this process is the unknown aliasing order, n. The traditional method of fluctu- 
ation spectrum analysis, the longitude-resolved fluctuation (LRFS), is insensitive to the sign 
or direction of drift, thereby doubling the number of aliasing possibilities in Eq. El DR99 used 
a new technique called the harmonic-resolved fluctuation spectrum, later shown to be equiv- 
alent to the two-dimensional fourier transform of the longitude-time dependence of intensity 
I Edwards fc Stappera , l2002i Edwards et al.l |2003) . Having eliminated the aliasing possibilities 



previously admitted by the LRFS, DR99 apparently conside red the aliasing question sol ved. 
However, even using the improved frequency measurements oflPeshpande fc R,ankinl l|2001[) . an 
inflnite range of possible values of n remain ( Edwards fc StappersL 120021) . 
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To obtain the other geometric parameters, IPeshpande fc RankirJ l|200l|) turned to polari- 
metric observations. Although the rotating; vec tor model has been shown to hold in only a 
minority of pulsars iJEverett fc Weisbera l200lj) . and B0943+10 exhibits a featureless linear 
sweep of position angle givin g no clear evidence for applicability of the model, it was taken by 
iDeshpande fc RankinI l)200l[l to apply. The values derived were reported to support the chosen 
aliasing solution, although as noted above other solutions are possible and only by combining the 
spect ral results and with polarimetry can a unique solution be obtained ( Edwards & Stapper^ 

Having obtained nominal values for all relevant geometrical parameters, 
IPeshpande & Rankin (2001) proceeded to form an animated map of the emission beam, 
a product which, if reliable, offers remarkable insights into conditions on the polar cap and/or 
in the magnetosphere. However, in my view this result is far from water tight. It relies on the 
reality of certain spectral features of marginal significance, the applicability of the rotating 
vector model for polarisation, and the applicability of the geometric predictions of the carousel 
model for drifting subpulses. Of the latter two conditions, the first has been proven to be 
generally false, and the second was until recently almost entirely untested. This was seen a 
sufficient motivation to rigorously test the geometrical predictions of the model. 

As I will outline below, the carousel model in its basic form has now been shown to be 
generally inapplicable. It is therefore advisable to approach any inferences made from pulsar 
data on this basis, including polar cap maps, with caution. 

3 TESTING THE CAROUSEL MODEL 

For reasons outlined in the preceding sections, a study was conducted of the subpulse phase 
variation in s everal pulsars. The firs t pulsar to be examined was PSR B0320+39, at a frequency 
of 328 MHz IjEdwards et al.L 120031 ). Surprisingly, the subpulse phase was shown to exhibit a 



phase "jump" of roughly 180°, near the centre of the profile (Fig.EJ. Such a feature cannot be 
reconciled with the predictions of the carousel model in its simplest form. Several features of 
the drift bands are suggestive of an origin in superposed, out-of-phase patterns: 

— the rapid 180° phase jump, 

— the vanishing subpulse amplitude at this point, 

— the continuity of the absolute value of the rate of change of subpulse amplitude at this 
point, and 

— the continuity of, and indeed, absence of any feature whatsoever in the total intensity profile 
at this point. 

Such features are familiar from other cases of superposed wave phenomena, for example in the 
presence of superposed orthogonal pol arisation modes of sh ifting dominance. 

Two possibilities were discussed bv Edwards et alJ l|2003() for the origin of the putative "dou- 
ble images" of the drifting subpulse pattern, both invoking an underlying circulating carousel 
system. The first suggests that the two images originate from either side of the magnetic pole, 
with that from the far side reaching the observer by virtue of magnetospheric refraction. For 
an odd number of beamlets and an axisymmetri c plasm a distr ibution, the images will be in 
antiphase. The second follows the suggestion of iRankiri Ijl99.'l) that emission occurs at two 
discrete heights in the magnetosphere. The divergence of field lines results in nested cones of 
emission, while aberration and retardation shift their centres somewhat. However, the differ- 
ence in emission height needed for sufficiently offset cones in this case is in excess of 10,000 km, 
which is strongly incompatible with other estimates of pulsar emission heights fe.sf. lPvks et al.l 
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Fig. 2 Subpulse phase jumps. (Left) Pulse profile (dotted line), subpulse amplitude 
(solid line) and phase (points, plotted twice with linear slope subtracted) for PSR 
B0320+39 at 328 MHz. (Right) Pulse profile (thick solid line), subpulse amplitudes 
and phases for orthogonal polarisation modes (thin dotted and dashed lines, slope 
subtracted), for PSR B0809+74 at 328 MHz (see ;Edwards 2004 for full description). 



The detection of the phase jump in PSR B0320+39 also led to the development of a com- 
pletely new model for drifting subpulses. Clemens fc Rosen (2004') suggested that drifting sub- 
pulses are a manifestation of non-radial oscillations of high wavenumber. In this model, the 
phase jump and amplitude nulling are neatly explained as corresponding to the passage of the 
line of sight over a nodal point in the oscillation. However, as a temporal oscillation, the subpulse 
phase should be a strictly linear function of pulse longitude, which is inc ompatible with minor 
but definite deviations from linearity detected by E dwards et al. [2003^. Although amplitude 
wind owing of the subpulses may produce apparently curved drift bands (Clemens &: Rosenl 
120041) . as noted in Sect. 1, the prime advantage of using subpulse phase to test models of 
drifting subpulses is that it is in fact unaffected by such windowing. 

Further studies revealed more complex deviations for the simple model predictions. 
lEdwards fc Stapoer j (J2003) found that at 1380 MHz, PSR B0320-h39 showed subpulses that 
were much weaker, with a phase slope that contained strange deviations from linearity and only 
weak evidence for a phase jump. Furthermore, observations of PSR B0809+74 revealed unex- 
plained deviations from the model curve at 328 MHz, which evolved to include a rapid phase 
jump of some ~ 120°at 1380 MHz, accompanied by an attenuation of amplitude. Analogous 
to similar polarimetric behaviour attributed to non-orthogonal polarisation modes, the latter 
was interpreted as suggestive of the superposition of subpulse patterns that were not in pure 
antiphase, and p otentially were not offset in p hase by a constant amount across whole the pulse 
longitude range. LEdwards fc StappersI l)2003(l suggested that such circumstances could also ex- 
plain the curious bidirectional subpulse drift seen in PSR B1918-I-19, in much the same way that 
polarisation position angle sweeps can undergo erratic jumps and sense changes in the presence 
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of no n-orthogonal radiation. Since then, further examples of phase steps llWeltevrede et all 



|2006|) and bidirectional drift l)Champion et all feoOSt IWeltevrede et all l200^ have been found. 

Observations made over thirty years ago have shown that the polarisation at a given 
pulse lo ngitude is periodi c ally m odulated between two orthogonal states, for at least two 
pulsars jManchester et al.lll975(l . Following a high resolution study of PSR 60809+74 by 
iR.amachandran et al. II2002I) . it was suggested that this effect could be explained via the super- 
position o f out-of-phase su bpulse patterns in the two orthogonal modes ( Edwar ds et all 1200.1 : 
iR.ankin fc Ramachandr anl l200,1l) . By decomposing the signal into its component orthogonal 
modes. Edwarda (;2004^1 showed that this was indeed the case in observations made at 328 MHz. 
Moreover, one of the modes showed a sharp phase jump accompanied by a reduction in subpulse 
amplitude (Fig. |21l, very similar to the feature seen in total intensity at 1380 MHz. The other 
mode showed no feature at all in this region, supporting the physical validity of the orthogonal 
mode decomposition. The implications of this result are startling: if the orthogonal modes are 
out of phase due to superposition then there are at least two images, but if the phase jump in 
one of the modes is also due to out-of-phase superposition, then there must be at least three 
distinct, out of phase "images" of the subpulse pattern, with polarisation effects somehow tied 
in to the multiple imaging process. 

Single pulse polarimetry of PSR B0809-I-74 at 1380 MHz, and observations of two other pul- 
sars at lower frequencies revealed significantly more complex behaviour. In these observations, 
the polarisation vector (composed of Stokes Q, U and V) was found to move along a periodic 
locus at the frequency of the subpulse modulation, however the first harmonic of this locus 
was elliptical, rather than linear as expected under the superposition of orthogonal polarisation 
modes. 

4 SUMMARY AND CONCLUSIONS 

Subpulse phase analysis presents a sensitive means for testing strong geometric predictions 
made by the standard carousel model for drifting subpulses. Testing this model is important 
because a variety of inferences have been made from observational results under the assumption 
that the model applies. In particular, procedure of "mapping" the polar cap of PSR B0943+10 
relies heavily on this model, as well as assuming, despite the odds, that the geometric model 
for polarisation applies to this pulsar. A number of studies of subpulse phase have now been 
conducted, with results that strongly disagree with the standard model. Several features of the 
observed phase slopes are suggestive of the superposition of out of phase "images" of subpulse 
patterns. A polarimetric study of PSR B0809-I-74 indicates that out of phase, orthogonally 
polarised subpulse patterns are observed in superposition, and that these patterns themselves 
can consist of the superposition of two or more out of phase components. At other frequencies, 
complex, non-orthogonal periodic modulations of the polarisation were seen, with unknown 
origin. 

In light of the strong deviations from the standard model observed in all cases so far, 
I suggest caution be exercised in the interpretation of observational results that assume its 
validity. 
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